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Abstract In our Paper I, we performed the spectroscopic observations of 1 1 confirmed 
GCs in M31 with the Xinglong 2.16m telescope and we mainly focus on the fits method 
and the metallicity gradient for the M31 GC sample. In this paper, we analyzed and dis- 

Q-i, cussed more about the dynamics, metallicity and age, and their distributions as well as 

the relationships between these parameters. In our work, eight more confirmed GCs in 

£h \ the halo of M3 1 were observed, most of which lack the spectroscopic information before. 

5/5 ■ These star clusters are located far from the galactic center at a projected radius of ~ 14 

to ~ 117 kpc, which are more spatially extended than that in the previous work. The Lick 
absorption-line indices and the radial velocities have been measured primarily. Then the 
ages, metallicities [Fc/H] and [a/Fe] have been fitted by comparing the observed spectral 
feature indices and the SSP model of Thomas et al. in the Cassisi and Padova stellar evo- 
lutionary tracks, respectively. Our results show that most of the star clusters of our sample 
are older than 10 Gyr except B290^ 5.5 Gyr, and most of them are metal-poor with the 
metallicity [Fc/H] < — 1, suggesting that these clusters were born at the early stage of 
the galaxy's formation. We find that the metallicity gradient for the outer halo clusters 
with r p > 25 kpc may not exist with a slope of —0.005 ± 0.005 dex kpc -1 and if the 
O 



outliers G001 and HI 1 are excluded, the slope dose not change significantly with a value 
of —0.002 ± 0.003 dex kpc -1 . We also find that the metallicity is not a function of age 
for the GCs with age < 7 Gyr while for the old GCs with age > 7 Gyr there seems to be a 
trend that the older ones have lower metallicity. Besides, We plot metallicity distributions 
with the largest sample of M3 1 GCs so far and it shows the bimodality is not significant 
and the number of the metal-poor and metal-rich groups becomes comparable. The spa- 
tial distributions shows that the metal-rich group is more centrally concentrated while the 
metal-poor group is occupy a more extended halo and the young population is centrally 
concentrated while the old populaiton is more extended spatially to the outer halo. 

Key words: galaxies: individual (M31) — galaxies: star clusters — globular clusters: 
general — star clusters: general 

1 INTRODUCTION 

One way to better understand the formation and evolution of the galaxies is through detailed studies of 
globular clusters (GCs), which are often considered to be the fossils of galactic formation and evolution 
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processes, since they formed at the early stages of their host galaxies' life cycles (iBarmby et al.Ll2000h . 
GCs are densely packed, very luminous, which usually contains several thousands to approximately 
one million stars. Therefore, they can be detected from great distances and are suitable as probes for 
studying the properties of extragalactic systems. Since the halo globular clusters (HGCs) are located far 
away from the galaxy center, they are very important and useful to study the dark matter distribution of 
the galaxy. Besides, since the HGCs are far from the galaxy center, the background of galaxy becomes 
much lower, which makes the observations much easier, compared to the disc GCs in the projected 
direction of galaxies. 

As the nearest (~ 780 kpc) and la rge spiral galaxy in our Lo cal Group, M31 (Andr omeda) contains 
a great number of GCs from 460 ± 70 dBarmbv & Huchrall200lb to ~530 (JPerinaet al.ll2010l) , which is 
an ideal laboratory for us to study the nature of the HGCs. A great many of new M3 1 HGCs have been 
discovered in t he recent yea r s, whi ch are important to study the formation history of M3 1 and its dark 



matter content. IHuxor et al. (|2004|) disco vered nine previously unknown HGCs of M31 using the INT 



survey. Subsequentlv. lHuxor et al.l (120051) found three new , extended GC s in the halo of M3 1 , which have 



characteristics between typical GCs and dwarf galaxies. iMackey et al.l (|2006J) reported four extended, 



low-surface-brightness clusters in the halo of M3 1 based on Hubble Space Telescope/Advanced Camera 
for Surveys (ACS) imaging. These star clusters are structurally very different from typical M3 1 GCs. 
On the other side, since they are old and metal-poor, they look like the typical Milky Way GCs. lHuxon 
(2007) found 40 new extended GCs in the halo of M31 (out to ~ 100 kpc from the galactic center) 
based on INT and CFHT imaging. These extended star clusters in the M31 halo are very similar to 
the diffuse star c lusters (DSCs) associated with early-type galaxies in the Virgo Cluster reported by 
iPeng et all d2006l) based on the ACS Vi rgo Cluster survey. In deed, the evidence shows that DSCs are 
usually fainter than typical GCs. Later, iMackev et al.l (120071) reported 10 outer-halo GCs in M31, at 
~15 kpc to 100 kpc from the galactic center, eight of which were newly discovered based on deep 
ACS imaging. The HGCs in their sample are very luminous, compac t with low metal licity, which are 
quite different from their counterparts in our Galaxy. More recently, iMa et al.l (12010b constrained the 
age, metallicity, reddening and distance modulus of B379, which also is an HGCs of M31, based on the 
multic olor photometry. 

In lFan et all J201ll) (hereafter Paper I) we observed 1 1 confirmed star clusters, most of which are 
located in the halo of M31, with the OMR spectrograph on 2.16m telescope at Xinglong site of National 
Astronomical Observatories, Chinese Academy of Sciences, in fall of 2010. We estimated the ages, 
metallicities, a-elements with the SSP models as well as the the radial velocities and they found that 
most of the halo clusters are old and metal-poor, which were supposed to be born at the early stage of 
the galaxy formation history. In this paper, we will continue the study of the HGCs of M3 1 with the 
same instruments and a larger sample. This allows us to be able to better understand the properties of 
the M31 outer halo. This paper is organized as follows. In S|2]we describe how we selected our sample 
of M31 GCs and their spatial distribution. In Sj3] we reported the spectroscopic observations with 2.16 
m telescope and how the data was reduced and the radial velocities and Lick indices were measured and 
calibrated. Subsequently, in [J4] we derive the ages, metallicities and a-element with x 2 — minimization 
fitting. We also discuss our final results on the metallicity distribution in the M31 halo. Finally, we 
summarized our work and give our conclusions in Sj5] 

2 SAMPLE SELECTION 

The sources were selected from the updated Revised Bologna Catalogue of M3 1 globula r clusters an d 
candidates (RBC v.4, available from htt p://www.bo.astro.it/M 3 1 jGalleti et"aTll2004ll2006ll2007l f2009). 
which is the latest and most comprehensive M31 GC catalogue so far. The catalogue contains 2045 
objects, including 663 confirmed star clusters, 604 cluster candidates, and 778 other objects that were 
previously thought to be GCs but lat er proved to b e stars, asterisms, galaxies, or Hll regions. Indeed, 
many of the halo clusters were from IMackev et all (120071) . who reported 10 GCs in the outer halo of 



M3 1 from their deep ACS images, of which eight were detected for the first time (see for details in [JT). 
In our work, our sample clusters are completely selected from RBC v.4. We selected the confirmed and 
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luminous clusters as well as being located as far as they could from the galaxy center, where the local 
background is too luminous for our observations. Finally, there are eight bright confirmed clusters in our 
sample, all of which are located in the halo of the galaxy. These clusters lack spectroscopic observational 
data, especially for the metallicity measurements. Thus it is necessary to observe the spectra of our 
sample clusters systematically and constrain the spectroscopic metallicities and ages in detail. 

The observational information of our sample GCs are listed in Table [T] which includes the names, 
coordinates, projected radii in kpc, exposures and observation dates. All the coordinates (R. A. and Dec. 
in Cols. 2 and 3) and projected radii from the galaxy center r p (Col. 4) a re all from RBC v.4 , which were 
calculated with M31 cente r coordinate 00 : 42 : 44.31 , +41 : 16 : 09.4 dPerrettet al.[l2002l) . PA = 38° 
and distance d = 785 kpc dMcConnachie et alll2005h . 

Table 1 The observations of our sample GCs. 



ID 



R.A. 

(J2000) 



Dec. 

(J2000) 



(kpc) 



Exposure 
(second) 



Date 



B289 


00:34:20.882 


+41:47:51.14 


22.65 


6000 


08/28/2011 


B290 


00:34:20.947 


+41:28:18.18 


21.69 


7200 


09/01/2011 


Hll 


00:37:28.028 


+44:11:26.41 


42.10 


5400 


09/01/2011 


H18 


00:43:36.030 


+44:58:59.30 


50.87 


5400 


08/29/2011 


SK108A 


00:47:14.240 


+40:38:12.30 


14.47 


3600 


08/28/2011 


SK112A 


00:48:15.870 


+41:23:31.20 


14.28 


5400 


08/29/2011 


MGC1 


00:50:42.459 


+32:54:58.78 


117.05 


3600 


08/28/2011 


H25 


00:59:34.560 


+44:05:39.10 


57.35 


5400 


09/01/2011 



We show the spatial distribution of our sample eight halo GCs and all the c onfirmed GCs from RBC 
v.4 in FigureQ] The large ellipse is the M31 disk/halo boundary as defined by lRacina(ll99lh . Note that 
all of our sample are located in the halo of M3 1 , which can help us to access the nature of galaxy halo 
with an enlarged cluster sample, compared to Fan et al.l (1201 lb . 



3 OBSERVATIONS AND DATA REDUCTION 

Our Low-resolution spectroscopic observations were all taken at the 2. 16m optical telescope at Xinglong 
Site, which belongs to National Astronomical Observatories, Chinese Academy of Sciences (NAOC), 
from August 28th to September 1st, 2011 (Please see Table [T). An OMR (Optomechanics Research 
Inc.) spectrograph and a PI 1340x400 CCD detector were used during this run with a dispersion of 200 
A mm -1 , 4.8 A pixel -1 , and a 3.0 "slit. The typical seeing there was ~ 2.5 ". The spectra cover the 
wavelength range of 3500 — 8100 A at 4 A resolution. All our spectra have S/N > 40. 

In order to calibrate our 2.16m data onto the Lick system, we also observed eleven Lick stan- 
dard stars (HR 6806, HR 6815, HR 7030, HR 7148, HR 7171, HR 7503, HR 7504, HR 7576, HR 
7977, HR 8020, HR 8165) near our field, which are selected from a catalogue of all 25 index measure- 
ments and coordinates for 460 stars (ref, available from http://astro.wsu.edu/worthey/html/system.html; 
IWorfhev & Ottavianilll997tlWorfhev et al.l[l994al) . Most of these standard stars are luminous (~ 5 — 6 
in V band), hence the exposure time we took was 20 second with the OMR system. 

The spectroscopic data were reduced following the standard procedures with NOAO Image 
Reduction and Analysis Facility (IRAF v.2.15) software package. First, the spectra have been bias and 
flat-field corrected, as well as cosmic-ray removed. Then the wavelength calibrations were performed 
based on Helium/ Argon lamps exposed at both the beginning and the end of the observations in each 
night. Flux ca librations were perfo rmed based on observations of at least two of the KPNO spectral 
standard stars dMassey et al.Ul988|) each night. The atmospheric extinction was corrected with the mean 
extinction coefficients measurements of Xinglong through the Beijing-Arizona-Taiwan-Connecticut 
(BATC) multicolor sky survey (H. J. Yan 1995, priv. comm.). 

Before the Lick indices were measured, the heliocentric radial velocities V r were measured by 
comparing the absorption lines of our spectra with the templates in various radial velocities. The typical 
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Fig. 1 Spatial distribution of M31 GCs. Our sample halo GCs are shown with filled circles 
and the confirmed GCs from R BC v.4 are mar ked with points. The large ellipse is the M31 
disk/halo boundary as defined in lRacind dl99ll) . 



internal velocity errors on a single measure is ~ 20 km s -1 . The estimated radial velocities V r with 
the associated uncertainties (Col. 2) are listed in Table|2] The published radial velocities V r (Col. 3) are 
also listed for comparisons. The systematic difference between our observed velocity and the catalogue 
velocity is found to be 29 ± 39 km s _1 and the standard deviation of the differences between our 
observed velocity and the catalogue velocity is 78 km s _1 for the five pairs of the radial velocities. 
It suggests that our measurements agree with those listed in RBC v.4 since the systematic difference 
between our measurements and the published values is not significant. 

Figure |2] shows the radial velocity V r (corrected for the systemic velocity of M31) as a function 
of the projected radii from the galaxy center. The Left panel is for the all confirmed GCs which have 
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Fig. 2 The distributions of radial velocity V r (corrected for the systemic velocity of M31). 
Left: all the confirmed GCs. Right: the HGCs only. 



the radial velocity V r measurements and the Right panel is for the HGCs, which refers to the GCs in 
the galaxy halo defined in Figure Q] It can be noted that the radial velocity distributions are basically 
symmetric in distributions either for all the confirmed GC sample or for the HGCs only. 

Table 2 The radial velocities V r of our sample GCs as well as 
the previous results. 



ID 


our work 


RBC v.4 


B289 


-96.81 ±47.27 


-181 ±30 


B290 


-488.73 ±43.14 


-381 ± 26 


Hll 


-173.02 ±39.63 




H18 


-300.48 ± 79.65 




SK108A 


-352.17 ±19.18 


-379 ± 38 


SK112A 


-342.68 ± 32.81 


-252 ± 46 


MGC1 


-412.67 ±17.13 


-355 ± 2 


H25 


-256.49 ± 55.28 





We plotted the radial velocities V r versus the projected radii r p in Figure|3]wher e the radial velocitie s 
have been corrected for the systemic velocity of M31 galaxy of 300 ± 4 km s~ 1 dPerrett et al.l [20021) . 
The Left panel is for all the confirmed clusters in RBC v.4 while the right panel is for the halo clusters 
which are defined in Figure Q] The points are the published measurements from RBC v.4 while the 
open triangles and the filled circles with errors are the measurements in Paper I and those in our work, 
respectively. In the Right panel, the symbols are the same as those in the Left panel. We find that the 
dispersion of the velocity becomes smaller when the GCs are locate further from the center of the galaxy 
with larger projected radius r v . It can be seen that the dispersion of the radial velocity becomes smaller 
when the projected radius r p is larger. 

Subsequently, all the spectra were shifted to the zero radial velocity and smoothed to the wave- 
length dependent Lick reso lution with a variable-width Gaussian kernel following the definition of 
IWorfhev & Ottavianil (1 1997b . i.e. 1 1.5 A at 4000 A, 9.2 A at 4400 A, 8.4 A at 4900 A, 8.4 A at 5400 
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Fig. 3 The radial velocity V r (corrected for the systemic velocity of M31) as a function of 
the projected radius. Left: all confirmed clusters and Right: the halo clusters. The filled circles 
with errors are the halo GCs from our sample while the points represent the velocities from 
RBC v.4 catalogue. 



A, 9.8 A at 6000 A. In deed, we measured all the 25 types of Lick indices stric tly by using the parame- 
ters and formulae from lWorthev et aLl(ll994al) and lWorfhev & Ottaviani d 19971). The u ncertainty of each 
index was estimated based on the analytic formulae (1 1)— (18) of ICardiel et al.l dl998l) . 

Eq.Q]is the linear fit formula for calibrating the raw measurements of our 2. 16m data to the standard 
Lick index system. The eleven standard stars are utilized for the fitting (Please see Figure |4]i and the 
results are listed in Table [3] 

EW ref = a + b ■ EW raw (1) 



4 FITTING, ANALYSIS AND RESULTS 
4.1 Model description 



iThomas et al.l J2003I) provided stellar population models including Lick absorption line indices for var- 
ious elemental-abundance ratios, covering ages from 1 to 15 Gyr a nd metallicities fr om 1/200 to 3.5 x 
solar abundance. These models ar e based on the stand ard models oflMaras ton ( 1998), with input stellar 
evolutionary tr acks fromlCassisi et al.l (1 19971) and [Bono et al J i ll 997b and a lSalpeterl (1 1955b stellar initial 
mass function. IThomas et al.l d2004b improved the models by including higher-order B aimer absorption- 
line indices. They found that these Balmer indices are very sensitive to changes in the a/ Fe ratio for su- 
perso lar metallicities. The lat est stellar popu l ation m ode l for Lick absorption- line indices (IThomas et al.l 
120101) is an improvement on IThomas et al.l (l2003b and IThomas et alj J2004J) . They were derived from 
the MILES stellar library, which provides a higher spectral resolution appropriate for MILES and SDSS 
spectroscopy, as well as flux calibration. The models cover ages from 0.1 to 15 Gyr, [Z/H] from —2.25 
to 0.67 dex, a nd [a/Fe] from —0.3 to 0.5 dex. In our work, we fitted our absorption indices based on 
the models o f Thomas et al.l d2010b . by using the two sets of stellar evolutionary tracks provided, i.e., 
ICassisi et all d 19971) and Padova. 
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Fig. 4 Calibrations of index measurements from the eleven standard stars of 2.16 m raw spec- 
tra with those from reference IWorthev & Qttavianil (1 1 997b : IWorthev et alj dl994al) . The linear 
fit coefficients of Eq. [3] have been derived to be used for calibrating our raw data to the Lick 
index system. 



4.2 Fits with stellar population models and the results 



Similar to lSharina et alj (120061) and our Paper I, the \ 2 —minimization routine was applied for fitting Lick 
indices with the SSP models to derive the physical parameters. As we measured 25 different types of 
Lick line indices listed in Table[3] all indices were used for the fitting procedure. As lThomas et al.lfcOlQl) 
provide only 20 ages, 6 metallicity [Z/H], and 4 a-element [a/Fe] for the SSP model, it is necessary to 
interpolate the original models to the higher-resolution models for our needs. We performed the cubic 
spline interpolations, using equal step lengths, to obtain a grid of 150 ages from 0.1 to 15 Gyr, 31 
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Table 3 The Linear Fit Coefficients a and b in Eq. [3] for trans- 
formations of the 2. 16m data to the Lick index system. 

Index a b 

H<5 A (A) -0.15 ±0.19 1.00 ±0.04 

HWp(A) 0.04 ±0.15 1.15 ±0.06 

CN1 (mag) 0.04 ± 0.01 0.84 ± 0.07 

CN2 (mag) 0.02 ± 0.01 0.98 ± 0.05 

Ca4227(A) -0.04 ±0.14 2.73 ±0.21 

G4300(A) -0.06 ±0.19 1.05 ±0.04 

H 7A (A) 1.73 ± 0.26 0.78 ± 0.03 

H 7F (A) 0.79 ±0.16 1.07 ±0.05 

Fe4383(A) -0.32 ± 0.36 1.46 ±0.10 

Ca4455 (A) 0.71 ± 0.56 1.50 ± 1.21 

Fe4531 (A) -0.30 ± 0.24 1.33 ± 0.09 

Fe4668(A) -0.16 ±0.31 1.16 ±0.06 

H/3(A) 0.17 ±0.16 1.03 ±0.05 

Fe5015 (A) -0.34 ± 1.07 1.44 ± 0.30 

Mgl(mag) 0.03 ±0.01 1.18 ±0.07 

Mg2 (mag) 0.03 ± 0.01 1.04 ± 0.03 

Mg&(A) -0.12 ±0.15 1.09 ±0.04 

Fe5270(A) -0.25 ±0.11 1.21 ±0.05 

Fe5335 (A) -0.04 ± 0.06 1.23 ± 0.03 

Fe5406(A) -0.10 ±0.08 1.39 ±0.07 

Fe5709(A) 0.11 ±0.03 1.30 ±0.06 

Fe5782(A) 0.12 ±0.10 1.41 ±0.24 

NaD(A) 0.21 ±0.36 0.91 ±0.14 

TiOl (mag) -0.07 ± 0.01 2.38 ± 0.19 

Ti02(mag) -0.07 ± 0.01 2.56 ±0.14 

[Z/H] values from —2.25 to 0.67 dex, and 5 1 [a/Fe] from —0 - 3 to 0.5 dex, which could make the fitted 
results smoother and more continuous. Since lWorthevl(ll994bb : lGalleti et alj (|2009) pointed out the age- 
metallicity degenaracy for most of the spectral feature indices measurements, which almost remain the 
same when the percentage change Aage/AZ = 3/2. Therefore, it is necessary for us to constrain the 
metallicity with the metal-sen s itive in dices before the fits. 

Fortunately, iGalleti et alj (120091) provide two ways to measure the metallicity from the metal- 
sensitive spectral indices directly. One method is through combining the absorption line indices Mg 
and Fe, [MgFe], which is defined as [MgFc] = ^Mgb (Fc), where (Fc) = (Fc5270 ± Fc5335)/2. 
Thus, the metallicity can be calculated from the formula below, 

[Fc/H] [MgFo] = -2.563 ± 1.1 19 [MgFe] -0.106 [MgFe] 2 ± 0.15. (2) 

The second way to obtain the metallicity from Mg2 is using a polynomial in the following, 

[Fe/H] Mg2 = -2.276 ± 13.053Mg2 - 16.462Mg2 2 ± 0.15. (3) 

Finally we obtained [Fe/H] avg with uncertainty in Table @] which is an average of the metallicities 
derived from the metallicity Eqs. |2]and|3] respectively. The averaged metallicity [Fe/H] avg will be 
used to const r ain the metallicity in the fits to break the age-metallicity trends/degeneracy. However, 
iThomas et al.l d2010h model only provide the metallicity parameters with [Z/H] and [a/Fe], thus we 
need to find a relationship between the iron abundance [Fc/H], total metallicity [Z/H] and a-element 
to iron ratio [a/Fe], which we can replace [Fe/H] with [Z/H] and [a/Fc] in the fit procedure. In fact, 
IThomas et al.l (120031) give the relation in Eq. [4] 

[Z/H] = [Fe/H] ± 0.94[a/Fe] (4) 



Name 


[Fe/H] avg 


B289 


-1.83 ±0.27 


B290 


-0.56 ± 0.63 


Hll 


-0.49 ±0.58 


H18 


-1.35 ±0.65 


SK108A 


-2.35 ± 0.22 


SK112A 


-1.62 ±0.43 


MGC1 


-2.06 ±0.33 


H25 


-2.74 ±0.47 
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Table 4 The metallicities [Fc/H] derived from the spectral in- 
dices [MgFe],Mg2. 



Notes: Here we define [Fc/H] avg = ' Fc / H ][M E F^+[Fe/H] Mg 2 

Here we would like to draw reader's attention that although the metallicity [Fe/H] has been deter- 
mined primarily, there are still many different ways to combine [Z/H] and [a/Fe] in the parameter grid 
of the model. Therefore, we still need to fit the age, [Z/H] and [a/Fe] simultaneously. Here, we would 
like to constrain the metallicity in the fits for |[Fc/H]fj t — [Fc/H] avg | < 0.3 dex, which is the typical 
metallicity uncertainty for the observations and it will make the fits more reasonable. Like the Paper 
I, the physical parameters ages, metallicities [Z/H], and [a/Fe] can be determined by comparing the 
interpolated stellar population models with the observational spectral feature indices by employing the 
X 2 —minimization method below, 




L%° del (age,[Z/H],[q/Fe]) \ 2 



(5) 



where L™ odol (age, [Z/H], [a/Fe]) is the i th Lick line index in the stellar population model for age, 
metallicity [Z/H], and [a/Fe], while L^ bs represents the observed calibrated Lick absorption-line in- 
dices from our measurements and the errors estimated in our fitting are given as follows, 

G i = °obs,i + ^modol.i- (6) 

Here, Cobs.i is the observat ional uncertainty while <7 mo dcU is the uncertainty associated with the mod- 
els of iThomas et alj J2010I) . These two types of uncertainties have been both considered in our fitting 
procedure. 

Ta ble [5] lists the fitted ages, [Z/H] and [a/Fe] with different evolutionary tracks of ICassisi et al.l 
( 1 19971) and Padova, respectively. In addition, we calculated the [Fc/H] cass ; s i and [Fc/H] pa d va by ap- 
plying the Eq.|4]to the fitted [Z/H] and [a/Fe]. For the reason of keeping consistency with Paper I, we 
adopted the metallicity [Fe/H] cass j s ; in the following statistics and analysis. From T able [5] we found 
that the ages, [Z/H] and the a-element [a/Fe] fitted from either ICassisi et all (1 1997b or Padova tracks 
are consistant with each other. Besides, it is worth noting that all of our sample halo GCs are older than 
10 Gyr in both evolutionary tracks except B290 (5 .5 to 5.8 Gyr), which is older than 2 Gyr and it should 
be identified as the "old" in ICaldwell et al.l J2009I) . Thus, it indicates that these halo clusters formed at 
the early stag e of the galaxy forma tion, which agrees well with the previous findings. 



Actuallv. lMackev et al] ( l2010l) conclude that the metal abundance of MGC1 is about [Fe/H] = -2.3 



and age is 12.5 to 12.7 Gyr through the color-magnitude diagram fitting. The age estimated agree well 
with our results while the metallicity is lower than ou r estimate [Fe/H] a y g = — 2.06 ± 0.33 in Table|4]or 
[Fe/H] cassis = -1.76±0.16 in Tabled Nevertheless jAlves-Brito et al.l(l2009l) found that the metallicity 



[Fe/H] = —1.37 ± 0.15 by combining the spectroscopic data and the photometric data, which is higher 
than our estimate. Hence, it can be seen that our result is just between the two results, suggesting that 
our result agrees with the previous conclusions. 

4.3 Metallicity Properties of Outer Halo 

The metallicity gradient of the halo star clusters and stars are important to the formation and enrichment 
processes of their host galaxy. Basically, there are two possible scenarios for the galaxy formation. 
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Table 5 


The % 2 — minimization 


Fitting Results Usin 


g 








iThomas et al 


J d2010t) Models with 


Cassisi et al 


J (119971) and 








Padova Stellar Evolutionary Tracks, respectively. 












Cas 


sisi 






Padova 




Name 


Age 


[Z/H] 


[a/Fc] 


[Fe/H] 


Age 


[Z/H] 


[a/Fe] 


[Fe/H] 




(Gyr) 


(dex) 


(dex) 


(dex) 


(Gyr) 


(dex) 


(dex) 


(dex) 


B289 


10.75 ±4.15 


-1.67 ±0.23 


0.34 ±0.16 


-2.09 ±0.27 


11.70 ±2.80 


-2.07 ±0.18 


-0.12 ±0.18 


-2.13 ±0.25 


B290 


5.80 ±2.40 


-0.99 ±0.05 


-0.26 ±0.05 


-0.85 ±0.07 


5.50 ±0.40 


-1.33 ±0.38 


-0.26 ±0.05 


-0.85 ±0.39 


Hll 


13.75 ± 1.25 


0.09 ±0.32 


0.08 ± 0.05 


-0.19 ±0.33 


13.60 ±0.20 


-0.10 ±0.24 


0.00 ± 0.06 


-0.21 ±0.24 


H18 


13.45 ± 1.45 


-0.47 ±0.37 


0.48 ± 0.02 


-1.07 ±0.37 


13.60 ±0.20 


-0.50 ±0.24 


0.48 ± 0.02 


-1.07 ±0.24 


SK108A 


13.60 ±0.30 


-1.53 ±0.18 


0.28 ±0.22 


-2.09 ±0.28 


13.55 ±0.45 


-1.48 ±0.23 


0.27 ± 0.24 


-2.09 ±0.32 


SK112A 


11.10 ±3.90 


-1.33 ±0.38 


0.25 ±0.25 


-1.35 ±0.45 


11.70 ±3.30 


-1.51 ±0.47 


0.10 ±0.40 


-1.42 ±0.61 


MGC1 


13.30 ±0.80 


-1.39 ±0.14 


0.42 ±0.08 


-1.76 ±0.16 


12.90 ± 1.30 


-1.39 ±0.14 


0.42 ±0.08 


-1.76 ±0.16 


H25 


13.60 ±0.30 


-1.98 ±0.20 


0.50 ±0.00 


-2.45 ±0.20 


13.50 ±0.50 


-2.03 ±0.05 


0.50 ±0.00 


-2.45 ±0.05 



One is that the halo stars and clusters should feature large-scale metallicity gradients if the enrichment 
timescale is shorter than the collapse time, which may be due to the galaxy formation as a consequence 
of a monolithic, dissipative, and rapid co llapse of a single massive, nearly spherical, spinning gas cloud 
dEggen et al.Lll962tlBarmbv et al.ll2000l) . The other one is a chaotic scheme for early galactic evolution, 
when the loosely bound pre-enriched fragments merge with the protogalaxy durin g a very long period o f 
time, in which case a more homogeneous metallicity distribution should develop dSearle & Zinnl[l978h . 
However, most galaxies are believed to have formed through a combination of these scenarios. 

Ivan den Berghl(u969l) : lHuchra et al.l(ll982l) showed that there is little or no evidence fo r a general ra- 
dial m e tallicity grad i ent for GCs within a rad ius of 50 arcmin. However, studies including lHuchra et al.l 
dl99ll) ; |Perrett et al](l2002l) ; lFan et all (12008) support the possible existence of a radial metallici ty gradi- 
ent for the metal -poor M3 1 GCs, although the slope is not very significant. iPerrett et al.l (120021) suggest 
that the gradients is —0.017 and —0 .015 dex arcmin -1 for the full sample and inner metal-poor clus- 
ters. More recently. Fan et al.l (120081) found that the slope is —0.006 and —0.007 dex arcmin -1 for the 
metal-poor subsample and whole sample while the slope approaches zero for the metal-rich subsam- 
ple. Nevertheless, all these studies are based on GCs that are located relatively close to the center of 
the galaxy, usually at projected radii of less than 100 arcmin. Recently. iHuxor et al.l d201 ll) investigated 
the me tallicity gradient for 1 5 halo GCs to r p = 117 kpc with the metallicity derived from the CMD 
fittings iMackev et alj J20061 120071 l2Q10t) and the authors found that the metallicity gradient becomes 
insignificant if one halo GC H14 is excluded in th eir Figure 6. We found that our result is consistent well 
with the previous findings of IHuxor et al.l (l201ll) . In Paper I, we found the slope of metallicity gradient 
is —0.018 ± 0.001 dex kpc -1 for the halo clusters sample extended to r p ~ 117 kpc from the galaxy 
center. Further, the slope turns to be —0.010 ±0.002 dex kpc -1 if only considering the clusters r p > 25 
kpc. 

Since we have spectroscopic observations of eight more halo confirmed clusters, it is interesting 
to check if the metallicity distribution/spatial gradient would change with an enlarged halo clusters 
sample. For the new observed data, as we recalled in iH.21 only MGC1 have the previous metallicity 
measurements from the literatures, which are very different for different works and our measurement is 
just the median value. Thus, we adopted our measurement. Finally we have a metallicity sample of 391 
entries in total. 

Figure [5] shows the metallicity as a function of projected radius from the galaxy center for all outer 
GCs with spectroscopic metallicity with r p > 25 kpc from the galaxy center. The slope of a linear fit 
is —0.005 ± 0.005 dex kpc -1 , which is marked with a solid black line. However, if the two highest 
metallicity star clusters G001 and Hll are excluded, the slope turns out to be —0.002 ± 0.003 dex 
kpc -1 , which is shown with the red dashed line. Thus, both of the cases suggest there is none metallicity 
gradient for the M3 1 outer halo clusters when r p > 25 kpc, which agree with the conclusion of Paper I. 
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Fig. 5 Metallicities [Fc/H] versus projected radii for the outer halo GCs with r p > 25 kpc 
from the center of the galaxy. The slope of the linear fitting is —0.005 ± 0.005 dex kpc -1 
(black solid line). However, if the two highest metallicity GCs G001 and HI 1 are excluded, 
the slope turns out to be —0.002 ± 0.003 dex kpc -1 (red dashed line). 



Table 6 The slopes of metallicity gradient by adding the random 
errors to the data. 



No. 


kail 


fc<-i 


1 


-0.013 ± 0.010 


-0.013 ±0.011 


2 


-0.003 ± 0.010 


0.000 ± 0.014 


3 


-0.011 ±0.012 


-0.008 ±0.012 


4 


-0.009 ± 0.009 


0.000 ± 0.012 


5 


-0.003 ± 0.013 


-0.036 ±0.021 


6 


-0.002 ±0.012 


0.004 ± 0.022 


7 


-0.004 ±0.013 


-0.002 ± 0.022 


8 


-0.009 ±0.012 


0.009 ± 0.015 


9 


-0.005 ±0.010 


-0.008 ±0.017 


10 


-0.013 ±0.011 


0.003 ± 0.008 



Therefore it seems that the "fragments merging" scenario dominated in the outer halo during the galaxy 
formation stage. 

It should be noted that the metallicity gradient is fitted based on the data of our observations and the 
literature and the metallicities from different literature m ay not the same. For instance, the metallicity of 
G00 l is [Fe/H] = -1.08 ± 0.09 in lHuchra et all (Il99ll) while [Fe/H] = -0.73 ± 0.15 in lGalleti et all 
(120091) . Thus we wonder how the slope would change when the data is changed. We simulated ten sets 
of random data from a = —0.5 to 0.5 and added them to the metallicities that we used in Figure[5]and 
then refit the slopes again for ten times separatedly and the results are shown in Table [6] It shows that 
the slope dose not change significantly when the simulated errors were added, suggesting that the slope 
is stable even the data from different measures. 
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Fig. 6 Metallicity [Fe/H] versus radial velocity V r (corrected for the systemic velocity of 
M3 1 ) for all the GCs with spectroscopic metallicities and radial velocity. The small points are 
from the literature; the squares are from Paper I; the triangles are from our measurement. 



Figure [6] shows the relationship between the metallicities and the radial velocities V r which 
have been correcte d for the systemic velocity of the M3 1 ga l axy. The spectroscopic metallicit i es are 
from the litera t ure dHuchra et ali Il99lt iBarmbv et all bOOOt iPerrett et all 120021 bailed et all 120091: 
ICaldwell et al.L 1201 ll) . Paper I as well as this work and the radial velocities V r are from the RBC v.4, 
Paper I and this work. It seems that there is no any relationship between the metallicities versus the 
radial velocities V r . 

Figure [7] shows the metallicities versus ages of the GCs. The metall i cities are from the literatur e 



dHuchra et al.Lll99ll:lBarmbv et alll2000HPerrett et al.Ll2002HGalleti et alll2009tlc"aldwell et alllJOllh . 
Paper I as well as this work and the ages are from the iFan et all (l2010h . Paper I and this work. We 
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Fig. 7 Metallicity [Fc/H] versus ages for all the clusters with spectroscopic metallicity and 
age estimates. The open triangles are the data from the literature; the filled circles are the data 
from Paper I; the filled triangles are the data from this work. The solid line represents the 
linear fit of GCs younger than 7 Gyr while the dashed line is the fit for the GCs older than 7 
Gyr. 



would like to see if there is any relationship between the ages and metallicities for these GCs. Actually 
we find that the relationships are different for the GC populations with different age. The slope of 
the GCs younger than 7 Gyr is k = 0.035 ± 0.021 while the slope of the GCs older than 7 Gyr is 
k = —0.095 ± 0.034, which is ~ 3er significant level. It suggests that for the GCs younger than 7 Gyr, 
there is no relationship between the age and metallicity while for the clusters older than 7 Gyr, it seems 
that the older GCs are more metal-poor (lower metallicity) and the younger GCs are more metal-rich 
(higher metallicity). 

Previously, many astronomers found the significa nt bimodal case in th e metallicity of M3 1 GC 
distribution by apply i ng the mixtu r e-mod el KMM test JAshman et all Il994l) . lAshman & Bird! dl993l) : 
iBarmbv et al.l (|2000|) ; iPerrett et al.l (|2002|) found the proportion of the metal-poor and metal-rich grou 
is ~ 2 : 1 to ~ 3 : 1 with the peak positions of [Fc/H] « — 1.5 and —0.6, respectively. iFan et al.l (120081 
examined the bimodality of metallicity distribution with a larger sample and the authors found the pro- 
portion is ~ 1.5 : 1 and the the peak pos itions are [Fc/H] ~ — 1.7 and ~ —0.7, respectively. However, 
the recent work of ICaldwell et al.l (1201 lb suggests that there is no significant bimodality or trimodality 
for metallicity distribution with a sample of 322 M3 1 GCs, most of which have spectroscopic metallic- 
ity with high S/N ratio. Since we have new observation data and a larger spectroscopic data sample, we 
are able to reexamine the bimodality of the metallicity distributions of M31 GCs. Figure [8] shows the 
metallicity distributions of the GCs and the HGCs, respectively. In the Left panel , the sample includes 
all the GCs which have spectroscopic metallic i ty from the literature dHuchra et al.l Il99ll IBarmbv et al.l 
l2000t IPerrett et al.l 120021 iGalleti et al.l 120091 ICaldwell et all 1201 lb and Paper I as well as this work. 
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Fig. 8 Metallicity distributions with bin size of 0.3 dex. Left: all the GCs with spectroscopic 
metallicities. The mixture-model KMM test was applied to divide them to two groups. Right: 
all the HGCs with spectroscopic metallicities. 



In total, there are 386 GCs with spectroscopic metallicity in the distribution. We applied the mixture- 
model KMM algorithm to the dataset and it returns an insignificant bimodality with p-value = 0.369, 
which means that a bimodal distribution is preferred over a unimodal one at 63.1% confidence level. 
The numbers of the metal-poor group and the metal-rich group are iVl = 196, N2 = 190, respectively 
and the mean values of the two groups are [Fe/H]i = -1.43 (a\ = 0.327) and [Fe/H] 2 = -0.73 
{<j\ = 0.215), respectively. As we can see from the plot, the proportion of the metal-poor and metal-rich 
group is ~ 1 : 1, which is lower than the published results. The reason why the bimodal case becomes 
more insignificant with larger sample may be that more intermediate metallicity GCs, which is between 
the two metallicity peaks, have been discovered and those intermediate metallicity GCs cause the dis- 
tribution to be unlikely a bimodal or trimodal distribution. Therefore, the previous works found that the 
metallicity distributions of M3 1 GCs is like that of the Milky Way and more recent works with more 
data show that they are less similar to each other, which may indicate that the formations of the two GC 
system was substantially different. In the Right panel, it show the metallicity distribution of the HGCs 
and obviously the metal-poor GCs dominate in the distribution. 

As the M3 1 GCs have been divided into two different groups by the KMM test in the metallicity dis- 
tribution of Figure [8] we would like to examine the spatial distributions of the two groups with different 
metallicity. Figure |9]plots the spatial distributions of the metal-rich and metal-poor groups. Note that the 
metal-poor group appear to occupy a more extended halo and much more widely spatially distributed 
while the metal-ric h group is more c entrally concentrated, which is consistent with the conclusions of 
IPerrett et al.l (l2002l) : lFan et al.l (120081) . 

Since we have the age estimates of the halo GCs in M31, we are curious about whether the spatial 
distributions of the young and old populations are the same or not. Here we used the de finition of "old 
popula tion" for age > 2 Gyr and the "young population" for age < 2 Gyr as that did in Caldwell et al.l 
d2009h . For the purpose of enlarging our sample, the age estimates for M31 GCs in lFan et all (1201 Oh 
and Paper I are also merged into our sample. Figure [10] plots the young and old population spatial 
distributions, respectively. It is obvious that the young population is more centrally concentrated and it 
traces the disk shape of the galaxy well. However, the spatial distribution of the old population is more 
dispersive and it seems that they do not trace the disk shape of the galaxy. 
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Fig. 9 The spatial distributions of HGCs with different metallicities. Left: metal-rich GCs; 
Right: metal-poor GCs. The two groups were divided by the KMM test of Figure|8] 
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Fig. 10 The spatial distributions of HGCs with young and old populations, respectively. Left : 
young clusters with age < 2 Gyr; Right: old clusters with age > 2 Gyr. 



5 SUMMARY AND CONCLUSIONS 



This is the second paper of our serial works for M31 halo globular clusters. In Paper I, we mainly 
focus on the fits method and the metallicity gradient for the M3 1 GC sample. In this paper, we focus 
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on the dynamics, metallicity and age, and their distributions as well as the relationships between these 
parameters. 

We selected eight more confirmed and bright GCs in the halo of M3 1 from RBC v. 4 and observed 
them with the OMR spectrograph on 2.16 m telescope at Xinglong site of NAOC in the fall of 2011. 
These star clusters are located in the halo of galaxy at a projected radius of ~ 14 to ~ 117 kpc from the 
galactic center, where the sky background is dark so that they can be observed in high signal-to-noise 
ratio. 

For all our sample clusters, we measured all 25 L ick absorption-line indices (see the definitions 
in, IWorthev et"aIT, Il994at IWorthev & Ottavianil 1 1997b and fitted the radial velocities. We found that 
distributions of the confirmed GCs and the halo GCs are basically symmetric to the systematic velocity 
of the galalxy. 

Similar to ISharina et al.l (120061) and our Pap er I, we applied the % 2 —minimization method to fit 
the Lick absorption line indices with the updated iThomas et al.l (1201 Ol) stellar population model in two 
stellar evolutionary tracks of Cassisi and Padova, separately. The fitting results show that most of our 
sample clusters are older than 10 Gyr except B290~ 5.5 Gyr and most of them are metal-poor with 
metallicity [Fc/H] < — 1 dex except HI 1 and HI 8, suggesting that these halo star clusters were born at 
the early stage of the galaxy's formation 

Again, we would like to study the metallicity gradient of the halo GCs by merging more spectro- 
scopic metallicity from our work, Paper I and the literature . We only considered outer halo clusters with 
r p > 25 kpc and the fitted slope is —0.005 ± 0.005 dex kpc" 1 . However, if two metal-rich outlier 
clusters G001 and Hll are excluded, the slope is —0.002 ± 0.003 dex kpc -1 , which does not change 
significantly. Furthermore, in order to eliminate the effect the errors of different observations, we added 
the random errors from a = —0.5 to 0.5 to the data and refit the slope agian for ten times. The result 
shows that the simulated errors do not affact the slope much. Thus it seems that met allicity gradient for 
M3 1 outer halo clusters dose not exist, which agrees well with the previous findings dHuxor et al.l 1201 lb 
and Paper I. This result may imply that the "fragments merging" scenario is dominated in the outer halo 
of the galaxy beyond 25 kpc from the center during the early stage of the galaxy formation. 

We do not find a relationship between metallicity and the radial velocity for M31 GCs sample. It 
seems that the metallicity is not a function of age for the GCs with age < 7 Gyr while for the old GCs 
with age > 7 Gyr t here seems to be a trend that the older ones have lower metallicity. This conclusion 
is similar to that of lFan et aU J2006b . who found a possible general trend of the age-metallicity relation 
with a large scatter. In addition, we plot metallicity distributions with the largest sample of M3 1 GCs so 
far and it shows the b imodality is not significant compared to the previous work. This is also found by 
IC aid well et al.l (l201lb . who used the newly observed spectroscopic data. We also find that the number of 
the metal-poor and metal-rich groups becomes comparable while the previous works show that the num- 
ber of metal-poor group is more than that of t he metal-rich one. This may be due to many intermedate 
metallicty metallicity of lCaldwell et al.l J201lb have been merged into our sample for our statistics. The 
spatial distributions shows that the metal-rich group is more centrally concentrated while the metal-poor 
group is occupy a more extended halo and the young population is centrally concentrated while the old 
populaiton is more extended spatially to the outer halo. This is easy to be understood as the old GCs are 
usually metal-poor especially for the halo GCs of M3 1 . 
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